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SECTION  I 


INTRODUCTION 

Parallel-plate  transmission  lines  are  often  employed  as  a guiding 
structure  for  electromagnetic  pulse  (EMP)  simulators.  This  type  of  transmission 
line,  as  shown  in  figure  1,  has  in  addition  to  the  dominant  TEM  mode  an  infinite 
number  of  higher-order  modes.  Because  of  the  open  nature  of  the  structure,  the 
propagation  constants  of  these  higher-order  modes  are  usually  complex, 
representing  the  propagation  as  well  as  radiation  loss  of  the  modes.  When 
an  object  to  be  tested  is  placed  in  the  present  structure,  it  is  illuminated 
by  electromagnetic  waves  consisting  of  a combination  of  TEM  and  a number  of 
higher-order  modes.  The  scattered  field  is  also  a superposition  of  TEM  and 
higher-order  modes.  Hence,  the  nature  of  the  higher-order  modes  in  the 
parallel-plate  transmission  line  is  worth  investigating  and  the  develop- 
ment of  a computer  program  to  obtain  the  propagation  constant  and  the  modal 
field  distribution  is  important  for  EMP  studies. 

The  problem  of  parallel-plate  transmission  lines  has  beer,  studied  by 
several  workers  (refs.  1 and  2)  using  various  approaches.  In  this  paper, 
a new  method  is  presented  for  attacking  this  problem.  The  method  is  an 
extension  of  Galerkin's  procedure  in  the  spectral  domain.  The  original 
version  of  this  method,  which  has  been  applied  to  many  mlcrostrip-type 
transmission  line  structures  (refs.  3 and  4),  is  extended  here  to  apply  to 
the  structure  with  complex  transverse  propagation  constants.  In  the 
following  sections,  the  formulation  of  the  problem,  numerical  procedures, 
and  some  results  are  presented. 
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Figure  1.  (a)  Cross  Section  of  a Finite-Width  Parallel-Plate  Waveguide 

(b),  (c)  Equivalent  Structure. 


SECTION  II 

FORMULATION  OF  THE  PROBLEM 

The  cross  section  of  the  parallel-plate  transmission  line  is  shown  in 
figure  1.  Assume  that  the  plates  are  infinitely  thin  and  that  both  the  plates 
and  the  medium  are  lossless.  It  is  well  known  that  the  dominant  mode  in  this 
structure  is  TEM  with  the  propagation  constant  identical  to  the  free-space 
wave  number.  All  of  the  higher-order  modes  may  be  classified  into  the  two 
sets,  TM  and  TE,  with  respect  to  the  z-direction.  Assuming  exp(jut  - j3z) 
variation,  where  the  propagation  constant  3 is  complex,  in  general,  all  of 
the  field  components  of  the  TM  and  TE  modes  may  generally  be  expressed  in 
terms  of  scalar  potentials  as 


where  k = 2ti/A  is  the  free-space  wave  number,  and  and  Pq  are  the  permittivity 
and  permeability  of  free  space,  respectively.  The  common  factor  exp(jo)t  - j3z) 
has  been  and  will  be  omitted  throughout  this  paper. 


7 


In  addition  to  the  TE  and  TM  classification,  the  symmetry  of  the  structure 


I 


i 

! 


allows  us  to  further  subclassify  the  modal  spectrum.  For  Instance,  the 
symmetricity  of  the  direction  of  the  z-component  of  induced  current  on  the 
plates  allows  us  to  subclassify  the  spectrum  into  the  four  cases  listed  in 
table  1.  The  detailed  formulation  will  now  be  given  for  Cases  1 and  3 of 
TM  modes  only  (odd  TM  modes) , and  only  the  resultant  equations  will  be 
summarized  for  the  rest  of  the  cases.  Because  of  the  nature  of  the  present 
method  of  analysis,  the  distinction  between  Cases  1 and  3 (and  between 
Cases  2 and  4)  is  not  necessary  in  the  formulation  process. 

Table  1 

CLASSIFICATION  BY  SYMMETRICITY 

Case  x>0,y>0  x<0,y>0  x<0,y<0  x>0,y<0 

1 + + 

2 + + + + 

3 + - - + 

4 + - + 

+ flows  in  the  positive  z direction 

flows  in  the  negative  z direction 

Such  distinction  is  undertaken  only  at  the  stage  of  preparation  for  numerical 
computation. 

In  the  odd  TM  mode  cases,  it  is  only  necessary  to  consider  the  equivalent 
structure  shown  in  figure  1(b)  where  the  y “ 0 plane  is  an  electric  conductor. 
Since  the  structure  is  infinite  in  the  x-direction,  the  electromagnetic 
boundary  value  problem  is  formulated  in  the  spectral  or  Fourier  transform 
domain  as  opposed  to  the  conventional  space  domain  formulation  (refs.  3 and  4). 


To  this  end,  let  us  define  the  Fourier  transform  $(a,y)  of  the  scalar  potential 
<|)(x,y)  via 


00 

4>j^(a.y)  = / 4>^(x,y)  e^“*  dx 

•»oo 


(3) 


where  i = 1 and  2 designates  the  regions  0 < y < d and  y > d,  respectively. 
The  transforms  of  field  components  may  be  defined  from  equations  (1)  and  (3) 
as 


>1 


= -ja  i^(a.y)  > 


“^0  3 r . . 


(4) 


I 


i 

t 


Since  <f>^  satisfies  the  wave  equation,  is  a solution  of 


where 


2 2 
Y = a 


^^(a,y)  = 0 


(5) 

(6) 


Because  of  the  boundary  conditions  = 0 at  y “ 0 and  the  radiation 

condition  at  y , the  solution  of  equation  (5)  is 


4'2^(a.y) 

= A (a) 

sinh  yy 

0 < y < d 

(7a) 

■ B(a) 

exp[-y(y  - d) ] 

y > d 

(7b) 

where  A and  B are  unknowns.  Note  that  Re  y > 0 and  Im  y > 0 are  to  be 
satisfied  so  that  equation  (7b)  represents  a valid  form  for  y -♦•  +»  . 
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The  next  step  is  to  apply  the  interface  conditions  at  y ■ d in  the 


transform  domain.  Since 

all  X 

E,l(x.d-) 

all  X 

Hxi(x,d-) 

lx|  > w 

the  interface  conditions  expressed  in  the  transform  domain  are 

E^^(a.d-)  = (8) 

ExiCa.d-)  = \2^a,d+)  (9) 

Hxi(a.d-)  - 


where  J is  the  transform  of  the  z-directed,  unknown,  induced  current  on  the 
z 

plate  at  y = d.  Substitution  of  equations  (4)  and  (7)  into  equation  (8)  gives 
the  relation  be:ween  A and  B.  Equation  (9)  is  automatically  satisfied  for 
A and  B.  If  these  quantities  are  substituted  in  equation  (10),  A or  B is 
expressed  in  terms  of  another  unknown  J^(a). 


where 


B(a)  = GqCu)  5^^“^ 
Gg(oi,B)  + coth  yd] 


(11) 

(12) 


Now,  the  final  boundary  condition  E^(x,:d)  = 0 for  |x|  < w is  imposed  in 
the  transform  domain.  Since  E^(x,d)  is  unknown  but  nonzero  for  |x|  > w,  it 
can  be  written  as 


10 


0 


E (x,d)  =/ 


2 2 

3 5 iiCx) 


B 


(13) 


Hence,  the  transform  is 


2 2 

E (a,d)  = j ^ — T-^  U(ot) 
z p 


(14) 


where 


-w 


U(a)  = / u(x)  e^°^*  dx  + / u(x)  dx 


w 


Elimixiating  B(ct)  from  equations  (8),  (11)  and  (14),  obtain 


GQ(a,B)  J^(a)  = U(a) 


(15) 


It  should  be  mentioned  that  equation  (15)  is  the  transform  of  the  integral 

equation  of  the  convolution  form  encountered  in  many  conventional  space  domain 

analyses.  It  may  also  be  worthwhile  to  mention  that  equation  (15)  contains 

two  unknowns,  and  U;  however,  it  is  possible  that  in  the  solution  one  of 

the  unknowns,  U,  may  be  eliminated  and  that  equation  (15)  is  solved  for 

J only, 
z 

Before  concluding  this  section,  let  us  summarize  the  resultant  equations 
for  other  symmetries  and  polarizations. 

TM,  even  in  y (Cases  2 and  4 in  table  1) 

G (a,B)  5 (a)  = U(a) 
e z 

G -g 

e (i)eQy[l  + tanh  yd] 


(16a) 

(16b) 
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TE,  odd  in  y 


d^(a)  = V(a) 


^0^“^  " ,„2  2 


JtL 


j(k  - g ) [1  + coth  yd] 


TE,  even  in  y 


L (a,g)  J (a)  = V(a) 
e z 


■ -;,2  ■„2 


j (k  - 3 )[1  + tanh  yd] 


For  TE  cases, 


-w 


V(a)  = / v(x)  dx  + / v(x)  dx 


w 


v(x)  = ^ H^(x,d) 


(17a) 

(17b) 


(18a) 

(18b) 


(19a) 

(19b) 
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SECTION  III 
METHOD  OF  SOLUTION 


In  this  section,  a method  of  solving  algebraic  equation  (15)  is 
discussed.  The  method,  which  is  applicable  to  solving  equations  (16a),  (17a), 
and  (18a)  as  well,  is  based  on  Galerkin's  procedure  applied  in  the  Fourier 
transform  domain. 

The  first  step  expands  the  unknown  J (a)  in  terms  of  known  basis 

z 

functions  J (a),  n = 1,  . . . ,N. 

n ’ * 

N 

J (a)  = y c J (a)  (20) 

z nn 

n=l 


where  c 's  are  unknown  coefficients  to  be  determined.  The  choice  of  J (a)'s 
n n 

is  such  that  they  are  the  Fourier  transforms  of  appropriate  functions  with 
finite  support,  viz.,  J^(x)'s,  the  inverse  transforms  of  J^(a)'s,  are  zero 
for  |x|  > w. 

Substituting  equation  (20)  into  equation  (15)  and  taking  an  inner 

product  of  the  resultant  equation  with  one  of  J^(a)'s,  m = 1,2 N, 

one  obtains 


where 


n=l 


mn 


c = 0 
n 


m = 1,2, 


.,N 


00 

K (S)  = / J„(u)  G (a,S)  J (a)  da 
mn  ' m O n 

«00 


(21) 


(22) 


The  right-hand  side  of  equation  (21)  is  zero  using  Parseval's  relation 


/ J (a)  U(a)  da 


m 


j(k‘  - 3‘) 


/ J (-x)  E (x,d)  dx  = 0 


because  and  are  nonzero  only  over  complementary  regions  of  x. 


% 
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When  J (x)'s,  whose  transforms  are  to  be  used  as  basis  functions  J (a)'s, 
n n 

are  chosen,  they  must  satisfy  certain  symmetry  requirements  in  addition  to 
being  zero  for  |x|  > w.  For  instance,  in  Case  1 in  table  1,  J^(x)'s  must  be 
symmetric  with  respect  to  the  y axis,  while  in  Case  3 they  are  required  to  be 
antisymmetric.  Furthermore,  it  is  desirable  to  use  J^(x)'s  which  well 
represent  the  edge  condition  at  x = ±w  where  the  actual  z-directed  current 
shows  square  integrable  singularity.  Using  these  basis  functions,  equation  (21) 
is  solved  for  unknown  propagation  constant  3 which  is  usually  a complex  number. 

After  3 is  obtained,  the  field  distribution  may  be  calculated  as  follows. 
Except  for  the  normalization  factor,  the  ratio  of  all  is  determined,  which 

gives  the  current  distribution 


J^(x) 


N 


I 

n=l 


c J (x) 
n n 


(23) 


The  field  distribution  of  may  be  obtained  from  equations  (7),  (11),  (15) 
and  (20) 


r 


E^(x,y)  = 


i ! ■=  f G„(a)  3 M da 

3 n •'  sinh  vd  0 n 

n=l  -«>  ' 

0 < y < d 

j - — S I ^n  ^ exp[-y(y  - d)]  Gp(a)  J^(a)  e 


(24a) 


S 


da  (24b) 


n=l 


y > d 


All  the  higher-order  TE  and  TM  modes  in  the  present  structure  may  be 
designated  by  a set  of  numbers  (p,q)  for  each  symmetric  subgroup.  The  number 
p is  associated  with  the  field  variation  in  the  x direction  and  number  q for 
the  variation  in  the  y direction.  In  the  present  method,  p can  be  predetermined 
by  the  appropriate  choice  of  basis  functions  for  the  current  distribution 


SECTION  IV 


NUMERICAL  PROCEDURE 


A numerical  algorithm  has  been  developed  for  TM  modes  of  Case  1 with 
the  mode  index  p = 0.  This  choice  of  p corresponds  to  the  higher-order 
modes  with  the  least  field  variation  in  the  x-direction.  The  basis  functions 
have  been  selected  so  that  the  qualitative  nature  of  the  actual  current  is 
well  represented.  Specifically  the  following  functions  have  been  employed: 


|x|  < w 
otherwise 


vO 


lx|  < w 
otherwise 


(25) 


(26) 


The  basis  functions  for  equation  (21)  are  the  Fourier  transforms  of 
equations  (25)  and  (26)>  which  are 


J2  (“) 


Jj^(a)  = 

4 sin  aw 
a 


2 sip  gw 
a 

1—  (1  - cos  aw) 

a w 


(27) 

(28) 


Although  any  number  of  functions  similar  to  these  given  by  the  above  equations 
could  be  used  as  a set  of  basis  functions,  N in  equation  (21)  was  set  to  two. 
Only  two  basis  functions  (27)  and  (28)  were  employed  for  economic  reasons, 
that  is,  for  minimizing  the  computer  time.  It  was  found,  however,  that  quite 
satisfactory  answers  were  obtained  for  many  microstrip  problems  by  using  only 
one  or  two  basis  functions  of  the  type  similar  to  equations  (27)  and  (28) 

(ref.  4). 


16 


Equations  (27)  and  (28)  are  substituted  into  equation  (22)  to  numerically 
compute  's,  and  complex  roots  S of 

det|K  (8)1=0  (29) 

are  found  by  a complex  root-seeking  algorithm.  This  algorithm  finds  the 
closest  zero  from  a given  starting  point.  At  this  stage,  only  the  mode 
number  p is  given  and  another  mode  number  q is  left  undecided.  It  is  possible 
to  correlate  the  value  of  q and  the  appropriately  chosen  starting  point  of  the 
algorithm. 

The  starting  values  of  the  route-seeking  routine  have  been  chosen  in  the 
following  way.  The  present  structure  can  be  viewed  as  a fictitious  closed 
waveguide  with  sidewalls  with  complex  surface  impedance.  For  well-guided 
modes,  these  sidewalls  may  be  very  close  to  the  magnetic  walls  since^ for 
such  modes,  the  radiation  loss  may  be  quite  small.  Hence,  the  propagation 
constant  8^  of  waveguides  with  magnetic  sidewalls  may  be  chosen  as  a starting 
point  of  the  algorithm.  (See  figure  2.)  For  the  present  mode  spectrum  of 
Case  1 and  TM,  8^  may  be  given  from  fl  ;ure  2a 


q = 1,  2,  . . . 

Note  that  we  are  concerned  with  p = 0 modes  in  the  present  numerical  computa- 
tions using  equations  (27)  and  (28)  for  basis  functions. 

It  is  hoped  in  the  numerical  algorithm  that  the  zero  of  equation  (29) 
closest  to  the  value  of  8^  be  obtained  for  a given  q;  such  a zero  is  called 


the  propagation  constant  8 of  the  TM  mode  of  Case  1. 

pq  pq 


Figure  2.  Classification  by  Symmetry  for  the  Determination 
of  the  Starting  Value  of  6. 
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SECTION  V 


NUMERICAL  RESULTS  ' 


Some  typical  computer  solutions  are  summarized  in  table  2 with  d,  w and 
the  wavelength  A being  input  parameters.  Also,  the  p = 0 and  q values  are 
specified. 

In  many  EMP  problems,  the  so-called  transverse  propagation  constant 
is  more  preferable  than  the  propagation  constant  3^^  (ref.  1).  The 
definition  of  a is 


pq 


(31) 


The  transverse  propagation  constant  here  corresponds  to  p^  given  in 
ref.  1.  Notice,  however,  that  unlike  in  ref.  1,  ct  as  well  as  S carries  two 
indices  p and  q because  the  structure  in  the  present  case  has  non-negligibie 
plate  width. 

The  magnitudes  of  the  E^  field  in  the  waveguide  cross  section,  which 
are  computed  using  equation  (24),  are  plotted  in  figure  3.  It  is  clear  that, 
the  field  decays  away  from  the  waveguide  in  both  the  x-  and  y-directions. 
Although  |e^|  must  be  zero  on  the  plate  at  y = d,  the  numerical  results  did 
not  predict  that  it  would  be  zero,  but  would  approach  zero.  It  is  hoped  that 
these  values  approach  zero  as  the  number  of  basis  functions  are  Increased. 

It  is  also  seen  that  the  number  of  peaks  in  the  y direction  for  0 < y < d 
coincides  with  the  given  value  of  q. 
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Table  2 


COMPUTER  SOLUTIONS  OF  PROPAGATION  CHARACTERISTICS 

(0,1)  mode 


w 

d 

d 

X 

2& 

; d 
s 

23  d 

pq 

2a  d 

pq 

No. 

Iter. 

2.5 

2 

24.32 

- j2.432 

21.76 

- 

j0.114 

-.1972 

+ 

jl2.58 

16 

10 

1 

10.88 

- jl.088 

10.82 

- 

j0.0652 

-.1104 

+ 

J6.382 

5 

5 

2 

24.32 

- j2.432 

21.72 

- 

j0.0448 

-.0728 

+ 

j 12. 644 

12 

3.33 

3 

37.17 

- J3.717 

20.82 

- 

j0.0335 

-.0222 

+ 

j 31.428 

45 

2.5 

4 

49.84 

- j4.984 

55.92 

- 

J6.68 

-27.28 

+ 

jl3.696 

15 

2 

5 

62.5 

- J6.25 

54.1 

- 

jO.163 

-.276 

+ 

j31.96 

24 

5 

10 

125.5 

- jl2.55 

107.2 

- 

jl4.72 

-22.28 

+ 

j70.80 

7 

5 

5 

62.5 

- j6.25 

59.6 

- 

jl4.5 

-28.46 

+ 

j30.36 

25 

5 

3.33 

41.4 

- j4.14 

45.8 

- 

jll.78 

-25.52 

+ 

j21.14 

29 

5 

2.5 

30.8 

- J3.08 

34.4 

- 

j3.94 

-15.92 

+ 

j8.52 

8 

5 

2 

24.4 

- J2.44 

21.8 

- 

j0.045 

-.078 

+ 

jl2.50 

12 

5 

1.67 

19.98 

- jl.998 

16.74 

- 

j0.060 

-.080 

+ 

jl2.58 

12 

5 

1.11 

12.46 

- jl.246 

12.36 

- 

j0.021 

-.040 

+ 

j6.50 

5 

5 

1 

10.88 

- jl.088 

10.78 

- 

j0.044 

-.074 

+ 

j6.46 

5 

2 

1 

10.88 

- jl.088 

10.29 

- 

j0.044 

-.063 

+ 

j7.22 

8 

1 

2 

24.32 

- j2.432 

19.28 

- 

j0.068 

-.081 

+ 

jl6.12 

15 

0.67 

3 

37.2 

- j3.72 

"k 

50.88 

(0,2) 

- jl2.66 
MODE 

-36.3 

+ 

jl7.70 

-7 

27 

10 

1 

0.2 

~ j0.02 

0.078 

jO.88 
X 10“ 

-7.6  X 10 

+ J9.56 

12 

*The  actual 

26  d was 

s 

0.0  - jO. 0, 

However,  23^d 

was  shifted  to  0.2 

- jO.C2 

for  numerical  convenience. 
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SECTION  VI 
CONCLUSIONS 


A simple  and  efficient  numerical  method  has  been  developed  for  analyzing 
a finite-width  parallel-plate  waveguide.  Sample  computations  based  on  this 
method  predicted  the  propagation  constants  of  the  modes  in  such  a structure 
and  the  field  plots  so  obtained  have  shown  the  expected  physical  nature  of 


these  modes. 
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